A gaussian model of the dynamics of an inextensible chain 
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In this work an approximated path integral model describing the dynamics of a 
inextensible chain is presented. To this purpose, the nonlinear constraints which 
enforce the property of inextensibility of the chain are relaxed and are just imposed 
in an average sense. This strategy which has been originally proposed for semi- 
flexible polymers in statistical mechanics, is complicated in the case of dynamics by 
the extra dependence on the time variable and by the presence of nontrivial boundary 
conditions. Despite these complications, the probability function of the chain, which 
measures the probability to pass to a given initial conformation to a final one, is 
computed exactly. The Lagrange multiplier imposing the relaxed condition satisfies 
a complicated nonlinear equation, which has been solved assuming that the chain is 
very long. 



INTRODUCTION 

The Rouse's [1] and Zimm's j^j models are able to describe satisfactorily the dynamical 
behavior of a polymer in a solution [3|]. There are however physical situations, such as for 
instance the case of a chain pulled by strong forces in experiments of micromanipulation of 
DNA {4-8], whose correct interpretation requires that the chain is inextensible. To imple- 
ment the property of inextensibility, it is necessary to introduce constraints in the stochastic 
equations that govern its dynamics. To this purpose, one may use generalized coordinates 
or Lagrange multipliers. Alternatively, it is possible to add forces, which become strongly 
repulsive whenever the chain attempts to abandon the region of the phase space in which the 
constraints are satisfied. Applying these methods one arrives at equations which determine 
the dynamics of a chain in the presence of constraints in a very rigorous way, but are also 
very complicated. Their solution, as much as the practical calculation of physical quanti- 
ties, require numerical simulations. An approach which leads to a simpler formulation of the 
problem consists in considering the chain as a system of particles held together by an elastic 
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potential. This potential must be chosen in such a way that in the equilibrium position the 
distance between two neighboring particles in the chain is equal to some constant a ^ 0. 
In the limit of infinite elastic constant, one obtains a discrete chain composed by particles 
connected together by massless segments of length a. In the limit in which a goes to zero, 
while the number of particles becomes infinite, one obtains the probability function of a 
continuous chain, which measures the probability that the chain passes after a fixed time 
tf from a given initial conformation to a given final conformation. It is possible to show 
that such probability function is equivalent to the partition function of a nonlinear field 
theory, which has been discussed in Ref. The appearance of nonlinearity is related to 
the inextensibility constraints which require that at any instant t and at each point of the 



chain, whose position is provided by the radius vector R(t, s), the relation 



8R(t,s) 
ds 



1 is 



the curve describing the spatial 
9j has a relatively simple formu- 



satisfied. Here we have denoted with s the arc-length o: 
conformation of the chain. The nonlinear model of Ref. 
lation in terms of path integrals and allows analytic calculations, such as for example that 
of the dynamical form factor in the semiclassical approximation 10]. Yet, it is important 



to have also a suitable approximation of that model which is able to simplify the functional 

— 1 . This delta function appears in the model as a result 



Dirac delta function 5 



dR(t,s) 



ds 



of the inextensibility constraints and it is a very nonlinear term. A similar delta function has 
been simplified so far in the case of the statistical mechanics of a freely hinged chain, where 
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;he time variable is not present, using the substitution S 



dR(s) 
ds 
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11] . In the latter formula, L represents the total length of the chain. Unfortunately, the 



extension of the approach of Ref. 11] to polymer dynamics is not feasible. To obtain a 
Gaussian approximation of the functional delta function comparable to that proposed by 
the authors of 11| also in the case of dynamics, we use here an alternative strategy bor- 
rowed from the statistical mechanics of semi- . The idea is to loose the 
condition that all arcs composing the chain must be inextensible and to require instead that 
the time averaged total length of the chain should fluctuate around a fixed average value L. 
This means that the length of the chain is allowed to change in time as it happens in the 
Rouse model, but only in such a way that its average is equal to L. This relaxed constraint 
is imposed in the probability function of the chain by introducing a Lagrange multiplier A, 
which is determined by minimizing the partition function as a function of A. With respect 
to the case of semi-flexible polymers, the situation is complicated by nontrivial boundary 
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conditions and the presence of the time. At the end we are able to compute exactly the 
explicit expression of the probability function, but for its minimization it is necessary to 
solve a complicated algebraic equation in A. We find its solution in the limit of very large 
values of the average total length of the chain L up to corrections of the order 4- included. 



THE GAUSSIAN APPROXIMATION 

We consider in this work the partition function of the model (9)]: 

pR(t f ,s)=R f (s) t . 

Z= VR{t, a )e- e f° dtS ° dsK2(t ' s) 5(i? /2 (t, s) - 1) (1) 

J R{0,s)=Ro(s) 

where R(t,s) : [0, t/] x [0,L] — > M. d is a two dimensional vector field with d components. 
The partial derivatives of R(t, s) with respect to t and s are denoted as follows: 

A( t , = = ««M ( 2 ) 

The parameter c appearing in the action (JI|) is given by: 

M 1 



(3) 



2L 2k B Tr 

ks denotes the Boltzmann constant, T is the fixed temperature of the thermal bath in which 
the chain fluctuates and r is the relaxation time which characterizes the rate with which 
the infinitesimal beads lose their speed due to friction. Finally, M is the total mass of the 
chain. 

It has been shown in jt], [3] that this model describes the dynamics of an inextensible 
continuous chain of length L obtained by performing the continuous limit of the partition 
function of a freely jointed chain. The dynamics of the chain is followed during the period 
of time t e [0,t/]. s G [0, L] is the arc-length measuring the distance along the chain. The 
vector field R(t, s) denotes the positions of the infinitesimal elements of length ds composing 
the chain at a given instant t. Z has the meaning of the probability function which measures 
the probability that the chain during its fluctuations passes from the initial conformation 
Ro{s) at the instant t = to the final conformation Rf(s) at the instant t = tf. The 
constraints enforcing the conditions that the length of the links connecting the beads should 
be constant become in the continuous limit the constraint: 

R' 2 (t,s) = l (4) 
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which is imposed in Eq. (JTJ) using a functional Dirac delta function. It is easy to check that 
this constraint implies that the length of every part of the chain is constant in time. Let us 
take for instance an arc of the curve R(t, s) delimited by the points R(t, s\) and R(t, sz). 
Its length €12 (t) is given by: in{t) = J*/ 2 ds\R'(t, s)\. Due to Eq. (jlj) it is possible to write: 
£i2{t) = ds — s% — s\. Thus 

iu(t) = (5) 

The probability function Z in Eq. jTJ should be completed by appropriated boundary 
conditions. We require that at the initial and final instants and tf the chain is in the 
configurations Ro(s) and Rf(s) respectively, i. e.: 

R{t f ,s) = R f (s) (6) 
R(0,s) = Ro(s) (7) 

Moreover, we suppose that the chain is closed, so that periodic boundary conditions will be 
chosen with respect to s: 

R(t, s) = R{t, s + L) (8) 

In the following, it will be convenient to expand the boundary conformations in Fourier 
series: 

+00 

R f (s) = Yl e2ninib n (9) 

n=— 00 
+00 

n=— 00 

a n and b n being constant vectors. The partition function of the model described by Eq. §T§ 
has been computed in joj in the semiclassical approximation. It is also possible to derive 
an expression of the dynamic structure factor of the chain always in the semiclassical ap- 
proximation as shown in [lo|. Despite these successes, the presence of the functional Dirac 
delta function in the right hand side of Eq. (CQ) makes the analytical treatment of the model 
complicated. For this reason, it would be nice to simplify it with a suitable approximation. 
For instance, in the case of statistical mechanics of polymers, where the time variable is not 



appearing, Edwards and Goodyear have shown for d = 3 that: 



5{R 2 {s)-a 2 )~e-ti ds ^ R ' 2 ^ (11) 
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While the extension of the results of 



llj to dynamics is not straightforward, it is however 



possible to relax the constraint (j3J) requiring that it is satisfied only in an average sense. To 
this purpose, following a procedure borrowed from the statistical mechanics of semi-flexible 

n nn 

polymers, see for instance [13] and [14J, [15] , we introduce a new parameter A and replace the 



probability function (JTJ with the following simplified one: 

~R(tf,s)=Rf(s) 
' R(0,s)=Ro(s) 

Next, the parameter A is determined by requiring that: 

dZ{\) 



r R(t f ,s)=R f {s) t . 

Z{\) = / VR(t,s)e-f° dtI ° ^ 2 M+A(K' 2 (t, s )-i)] (12) 

J R(0,s)=Rn(s) 



dX 

Clearly, the above equation is equivalent to: 

1 . 1 rL 



(13) 



dt- / dsR' 2 {t,s) ) = 1 (14) 
tf Jo L J I 

where 

r R(tf,s)=R f (s) t 
((...))= / VR ^ s)e -^dtf L ds [cRHt,s )+ X(R'Ht,s)^)] { _^ ) (15) 

J R(0,s)=R (s) 

In words, with (I14p the functional integration over the conformations R(t, s) has been ex- 
tended to curves whose arcs do no longer satisfy the constant length condition of Eq. (JHJ). 
Also the total length of the chain is allowed to change in time, but its average performed 
over all possible conformations and during the period of time tf must be equal to L. 



CALCULATION OF THE PROBABILITY FUNCTION Z(X) 

The rest of this work is dedicated to the derivation of the probability function Z(X) and 
of the parameter A, which up to now is the still unknown ingredient of Eq. ffl2|) . Unlike the 
case of semi-flexible polymers, we have to deal here with the nontrivial boundary conditions 
(J6j) and (JTj). To get rid of them, it is convenient to split the bond vectors R(t, s) into two 
components: 

R(t,s) = R d (t,s)+r(t,s) (16) 
R c i(t, s) is a solutions of the classical equations of motion related to the action: 



t f pL 

2u „\ , WD'2 



S = J dt J ds 
'o Jo 



cR 2 (t, s) + X(R ,2 -l) (17) 



6 



appearing in the probability function of Eq. (!T2|) : 

cR d + XR d = (18) 

R d (t,s) satisfies the boundary conditions ([6]) and (J7j): 

Rd(t f , s) = R f {s) R d (0, s) = R {s) (19) 

r(t, s) represents instead the fluctuations around the classical conformation R d (t, s). r(t, s) 
satisfies the Dirichlet boundary conditions: 

r(t f ,s) = r(0,s) = (20) 

With respect to the variable s both R d (t,s) and r(t,s) satisfy periodic boundary condi- 
tions analogous to those of Eq. (jSJ). After some calculations and without performing any 
approximation, it is possible to show that: 

Z{\) = e- s ^e Lt f x Z fluct (X) (21) 

where the action S(R d ) contains the contribution coming from the classical background 
conformations: 

S(R d ) = J f dtj ds(cR 2 cl + \rt£) (22) 
while Zfi uct (X) is the partition function of the fluctuations: 

Z fluct (X) = J Vre-ti fdt ti d < cr2+Xr ' 2 ) (23) 

First, we will compute S(R c i). It is possible to check after some integrations by parts that 

S{R cl )=cJ^ ds(R f {s)-R4t f ,s)-Ro{s)-R cl {0,sf) (24) 



with R cl (t f ,s) = and R d (0, s) = 

t=t f 



. The solution of the classical 

t=o 



equations of motion f JTB"]) satisfying the boundary conditions of Eq. ffl~9]) is: 
where 

a, = (») 
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Substituting the expression of the classical solution R d of Eq. ( )25|) in Eq. ( 12^|) . we obtain: 

1 



S{R C 



-cL P- 



sinh(/3_ n fy) 



i , v . cosh(/3_ n t / ) 

' °-nJ H -wo ; — T l°n ' °-n + a n ' a -nJ 



(27) 

Next, we derive the partition function of the fluctuations Zfi uct (X). From Eq. ( 123]) it turns 
out in d = 3 that: 

Z fluct (X) = (det'A)-f (28) 
where A is the differential operator: 



A-c— + A— 

5t 2 ds 2 



(29) 



The prime after the symbol of determinant in Eq. ( 1281) means that the zero modes of the 
operator A should be projected out from the expression of the determinant. At this point 
we compute the eigenvalues of A. To this purpose, we have to solve the equation: 



d 2 d 2 



(30) 



Dirichlet boundary conditions with respect to the time variable are assumed according to 
Eq. (|2Q|) . Moreover, ip should obey periodic boundary conditions with respect to s. It is 
possible to check that there exist solutions of Eq. ( |30|) with such boundary conditions when 
E takes the values: 



A c 
E n k = An' z I n 2 — + k 2 ^ 
1 L 2 -21} 



n,k = 0,±1,±2, 



(31) 



These solutions are the eigenfunctions: 

+oo 
n=—oo 

where 



ijj n ,k(f) = Asinh 



u 



4n 2 n 2 



A- 



L 2 



E 



n.k 



(32) 



(33) 



and A is a normalization constant. Excluding the trivial zero mode occurring for n = k = 0, 
we find that 



det'A = Yl' n . 



k>0 



(34) 
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In the above equation the prime in the product ]Tnfc>o means that the case n = k = 0, 
which corresponds to the excluded zero mode, should be omitted. Now we rewrite det'A as 
follows: 



A 



•I) J 2 J rin,/c>0 



An 2 n z + 



cL 2 k 2 



2t 2 X 



Using the (— function regularization it is possible to show that rin,fc>o"p- = 

(A) 3. Thus 

cL 2 k 2 



(35) 



det'A 



. 5 

A \ 3 
L 2 



n; 



n,k>0 



An 2 n 2 4 



2t 2 X 



(36) 



The remaining semi-infinite product has been evaluated by several authors, see for instance 



16N20|. Here we present just the result of the calculations 

cL 2 k 2 



n: 



n,k>0 



2t 2 X 



exp 



In r - — + In (1 
12 ^ V 



-1-KTlT ^ 



n>0 



where 



^ 2= 2gA 
cL 2 



(37) 



(38) 



Summarizing, the partition function of the fluctuations given in Eq. (128]) becomes 



Zfiuct(X) = exp 



15, t z c\ 3, 7TT 3v^, , 
— hi — y --IhtH > ln(l 



— 27rnr 



(39) 



// n>0 

Putting together the results of Eqs. ( 127|) and (|39|) . the total probability function Z(X) of 
Eq. ( !2~TT) may be rewritten as follows: 



Z{X) = e F ^ 



(40) 



where 



F(X) = -S(R d ) + Lt,A - ^ln(^)-^lnr+^-^ln(l - ) 



T 2 C 



TIT 3 



2t 



(41) 



n>0 



Clearly, the condition f lT3"|) that determines A is equivalent to the condition g£ = 0, i. e. 



dS(R cl ) 
dX 



- Ltt 



211 TTtf II 1 

~8~A ~ 8" 




L\ 2c y/x 



+ 2^~ 



2iin 



-1-KTlT 



n>0 




^ V 2 Cv ^ 



(42) 



As it stands, the above equation is too difficult to be solved analytically with respect to A. 
For this reason, we assume that the chain is very long, i.e. L » 1, so that it is possible 
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to expand A and the coefficients appearing in Eq. (142!) in powers of L 1 . In doing that, we 
suppose that A remains finite in the limit L — > +oo, so that: 

A ~ A + T Ai + . . . (43) 

This hypothesis on A, which will be checked a posteriori, has a physical motivation. In 
fact, if A becomes infinite with the average chain length L, then conformations for which 
jj f* f dtj^ Jq dsR' 2 < 1 will be strongly preferred in the probability function of Eq. ffl2|) . a 
fact which is clearly unphysical. 

At this point, we are ready to expand all the terms entering in Eq. (02]) in powers of j j . 
This expansion is almost trivial apart from the case of the two terms 

h = (44) 

and 

h = -T 2 ™ e- 2 ™ rtj -J^— (45) 

Ii depends on A implicitly through the coefficients (3 n , while I2 has an explicit dependence 
on A together with an implicit dependence through the coefficient r. Provided Eq. ( ]4"3"|) is 
fulfilled [211], we note that both parameters (3 n and r, being proportional to 4-, are going to 
zero when L goes to infinity. Using the above considerations and neglecting contributions 
of order jj or higher, we obtain: 

I x ~ (46) 

h - A^/cJ__A_i_ * [ItL (47) 
87rt / V2 A f 4A 167rt / V2 A f 

After completing the expansion in powers of \ of the remaining terms, Eq. fj42|) becomes: 

^1 A.Asi 9 , .... 




Sntf V 2 A | J J ' y 8 A l&irtf 




Equating to zero separately the coefficients accompanying different powers of L we obtain 
two relations which are able to determine both Aq and Ai: 



9c 

A, = A (50) 
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As a result it is possible to write down the expansion of A which is valid for large values of 
L up to the order j: 



This, together with Eqs. (140THTT) concludes the evaluation of the probability function (fl2l) . 



In this article a gaussian approximation of the nonlinear model ([1]) has been presented, 
which describes the probability of a chain to pass from a given conformation to another. 
The approximation consists in replacing the local constraint @ with a milder one, that 
requires only that during the fluctuations of the chain, its average length must be equal to 
L. To enforce this new constraint, the Lagrange multiplier A has been introduced. A has 
been determined exploiting the condition that, as a function of A, the probability function 
Z(A) of Eq. ( fl2l) has a minimum. Eq. ( fl4l) shows that this condition is equivalent to im- 
pose the relaxed constraint. While it is possible to compute the probability function Z(X) 
exactly despite the complications due to the presence of nontrivial boundary conditions, its 
minimization requires the solution with respect to A of the complicated algebraic equation 
given in Eq. fj42|) . This equation has been solved in the limit of very large values of L up to 
the second order included, see Eq. ( IBTT) . It is shown that A goes to zero for increasing values 
of the evolution time tf. This means that the deviations of the length of the chain from 
the expected value L increase with the increase of the time in which the chain is allowed to 
fluctuate. 

The present approach avoids the difficulties involved with the functional delta function of 
the original probability function (Cp) and it is still able to take into account the inextensibility 
of the chain, at least in an approximated way. The probability function Z(X) of Eq. ( fl2l) 
with the relaxed constraint has been computed here exactly, see Eqs. (HOTFfTl) . while the 
analogous probability function with the full constraint ([1]) could be derived only in semi- 
classical approximation. The approximated model of Eq. f fl2|) may be used as a basis for 
further computations of physical observables, such as for instance the dynamic form factor 
of a chain. 




(51) 



CONCLUSIONS 
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